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Euclidean dense matter generically suffers from the fermion sign problem. However, we argue that the sign 
problem is absent if one considers only low-energy degrees of freedom. Specifically, the low energy effective theory 
of dense QCD has positive Euclidean path integral measure, which allows one to establish rigorous inequalities 
showing that the color-flavor locked (CFL) phase is the true vacuum of three flavor, massless QCD. We then 
describe a method for simulating dense QCD on the lattice. We also discuss applications to electronic systems in 
condensed matter, such as generalized Hubbard models. 



Quark matter is described by a partition func- 
tion, 



dAdet (M) e 



-5(A) 



where [i is the chemical potential for the quark 
number and M = + is the Euclidean 

Dirac operator. Since the Dirac operator M con- 
sists of both anti-Hermitian and Hermitian oper- 
ators, it has in general complex eigenvalues. Fur- 
thermore, it is not related to its Hermitian con- 
jugate by any similarity transformations, 



M = 7 ££>£ 



(2) 



The determinant of M is therefore complex for 
generic gauge fields, which has thus far made lat- 
tice simulations very difficult |1I2I3) . 

However, it is shown here that the complex- 
ness of the measure of fermionic dense matter can 
be ascribed to modes far from the Fermi surface, 
which are irrelevant to dynamics at sufficiently 
high density in most cases, including quark mat- 
ter |4l5j . For modes near the Fermi surface, 
there is a discrete symmetry, relating particles 
and holes, which pairs the eigenvalues of the Dirac 
operator to make its determinant real and non- 
negative. 

As a simple example, let us consider fermionic 
matter in 1+1 dimensions, where non-relativistic 



fcrmions interact with a gauge field A. The action 
is in general given as 



(1) S= dTdx^[(-d T + t<f> + e F )-e{~id x +A)}^, (3) 



where e(p) ~ p 2 /(2m) + • • • is the energy as a 
function of momentum. Low energy modes have 
momentum near the Fermi points and have en- 
ergy, measured from the Fermi points, 

dE 

E[p±p F ) ~ ±v F p, v F = 



Op 



(4) 
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If the gauge fields have small amplitude and are 
slowly varying relative to scale p F , the fast modes 
are decoupled from low energy physics. The low 
energy effective theory involving quasi particles 
and gauge fields has a positive, semi-definite de- 
terminant. 

To construct the low energy effective theory 
of the fermionic system, we rewrite the fermion 
fields as 

i>{x,T) = ^ L {x,r)e +i ^ x + ^ R (x,T)e- l ^ x , (5) 

where ipL,R describes the small fluctuations of 
quasiparticles near the Fermi points. Using 

e ±t PF x E(-id x + A) e^P FX ip(x) w ± v F (-id x + 
A)ip(x), we obtain 

^eff = J dr dx [ip^ L (—d T + i(j> + id x — A)ip F 

+ tf R {-dr +%<!>- id* + A)$ R ]. (6) 
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Introducing the Euclidean (1+1) gamma matrices 
7o,i,2 andV'L.fl = ^(1±72)' ( A? we obtain a positive 
action: 



drdxip^id^+iA^ip^ / drdxipflip.{7) 



S c ff z 



Since ^) = 72-^72, the determinant of ^ is posi- 
tive, semi-definite. 

In this example, we see that low energy modes 
near the Fermi surface can be integrated (in favor 
of a determinant), leading to an effective theory 
without any sign problem whatsoever as long as 
they couple to slowly varying background fields. 
QCD at high baryon density falls into this cate- 
gory, since the coupling constant is small at high 
energy due to asymptotic freedom, and the corre- 
sponding fluctuations in the gauge fields at high 
density are small. 

A low energy effective theory of QCD at high 
density, called as High Density Effective The- 
ory (HDET), has been derived by one of us UJ. 
(Renormalization group analysis of the Fermi sur- 
face effective theory appears in [J].) Consider a 
quark in one of the patches that cover the Fermi 
surface only once. Its momentum can be decom- 
posed as 



ffy = A«V + In, |^| < A , Aj_ (< (j,), 



(8) 



where A and Aj_ are the sizes of patches perpen- 
dicular and parallel to the Fermi surface respec- 
tively, much smaller than the chemical potential 
but larger than the scale of interest. The normal- 
ization is enforced by a condition, 



V / d 2 l ± = 4tt P 



(9) 



patches * 

The modes near the Fermi surface are given as 



ip + (v F ,x) 



1 + a-VF ^-^gjr-x 



ip(x) 



(10) 



where a = 707 and Vf is the Fermi velocity of the 
modes. HEDT of quark matter is then described 
by 

Ajdet - i>+i^D^ + -^ +1 (l/) ± ) 2 ^ + + • • • , (11) 



where 7^ = (j°,vfVf '7*) = 7 M — 71 ■ 
that the leading term has a positive determinant, 



since 



Mn 



7l f • D(A) = 7 5 M e t ft75 



(12) 



In order to implement this HDET on lattice, it 
is convenient to introduce an operator formalism, 
where the velocity is realized as an operator, 



v = 



-i d_ 
V = V f d$ 



(13) 



Then the quasi quarks near the Fermi surface be- 
come 



, , , . v 1 + OL ■ V , 

tp = exp (+ifix ■ v) ip+. 



(14) 



Now, neglecting the higher order terms, the La- 
grangian becomes with X = cxp(i/ix ■ v)(l + a ■ 
«)/2, 



C 



HDET 



(15) 



where = X* X denotes soft gluons whose 
momentum \p^\ < \i. Since v ■ dv ■ 7 = d ■ 7 , we 
get 7y 1 (9 M = 7^9^, which shows that the operator 
formalism automatically covers modes near the 
full Fermi surface. 

Integrating out the fast modes, modes far from 
the Fermi surface and hard gluons, the QCD par- 
tition function Q becomes 



Z(ji) = J dA + det (M eff ) e" 



S of f(A + ) 



(16) 



where 

Soft = 

and A± 



1 M 2 

_pa pa ao, ao, 



(17) 



A — Au, the Debye mass M 



y / Nf/(2Tr 2 )g s n . At high density the higher order 
terms (~ A//i) are negligible and the effective ac- 
tion becomes positive, semi-definite. Therefore, 
though it has non-local operators, HDET in the 
operator formalism is free from the sign problem 
and can be used to simulate the Fermi surface 
physics like superconductivity. Furthermore, be- 
ing exactly positive at asymptotic density, HDET 
allows to establish rigorous inequalities relating 
bound state masses and forbidding the breaking 
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of vector symmetries, except baryon number, in 
dense QCD 0. 

With the help of the previous two examples, 
we propose a new way of simulating dense QCD, 
which evades the sign problem. Integrating out 
quarks far from the Fermi surface, which are sup- 
pressed by 1/fj, at high density, we can expand the 
determinant of Dirac operator at finite density, 



det (M) = [real, positive] 



1 



(18) 



As long as the gauge fields are slowly varying, 
compared to the chemical potential fj,, the sign 
problem can be evaded. As a solution to the sign 
problem, we propose to use two lattices with dif- 
ferent spacings, a finer lattice with a lattice spac- 
ing eidct ~ A* -1 f° r fermions and a coarser lattice 
with a lattice spacing a gauge -C /i -1 for gauge 
fields and then compute the determinant on such 
lattices. 

The determinant is a function of plaquettes 
{U x ^} which are obtained by interpolation from 
the plaquettes on the coarser lattice with spac- 
ing a gaU ge- To get the link variables for the finer 
lattice, we interpolate the link variables U XAI G 
SU(3) (see Fig.^J: Connect any two points 51,52 
on the group manifold as 



g(t)=gi+t(g 2 -gi) , 0<i< 1 



(19) 




Figure 1. Simulation with two lattices with dif- 
ferent lattice spacings 



For importance sampling in the lattice simula- 
tion, one can use the leading part of the deter- 
minant, [real, positive] . This proposal provides a 



nontrivial check on analytic results at asymptotic 
density and can be used to extrapolate to inter- 
mediate density. Furthermore, it can be applied 
to condensed matter systems like High-T c super- 
conductors, which in general suffers from a sign 
problem. 
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